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ABSTRACT 

We extend our previous analysis which used generalized luminosity functions (GLFs) 
to predict the number of quasars and galaxies in low-radio-frequency-selected samples 
as a function of r edshift, radio luminosity, narrow-e mission- line luminosity and type 
of unified scheme l)Grimes. Rawlings fc Willottll2004) . Our extended analysis incorpo- 
rates the observed submillimetre (850 /Ltm) flux densities of radio sources, employs a 
new method which allows us to deal with non detections, and focuses on the high- 
luminosity population. First, we conclude that the submillimetre luminosity Lgso of 
low-frequency-selected radio sources is correlated with the bolometric luminosity LboI 
of their quasar nuclei via an approximate scaling relation £350 oc L'^1\^^^'^- Second, we 
conclude that there is quantitative evidence for a receding-torus-like physical process 
for the high-luminosity population within a two-population unified scheme for radio 
sources; this evidence comes from the fact that radio quasars are brighter in both nar- 
row emission lines and submillimetre luminosity than radio galaxies matched in radio 
luminosity and redshift. Third, we note that the combination of a receding-torus-like 
scheme and the assumption that the observed submillimetre emission is dominated by 
quasar-heated dust yields a scaling relation Ls5o oc Lb^j which is within the errors of 
that determined here for radio-selected quasars, and consistent with that inferred for 
radio-quiet quasars by Willott. Rawlings & Grimes (2003). 
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1 INTRODUCTION 

Unified schemes for radio galaxies and radio quasars pro- 
pose that both types of objects are physically similar 
but viewed at different angles between the radio axis 
and the line of sight. A comprehensive introduction to 
th e subtleties involved in su c h sch emes can be found 
in iGrimes. Rawlings fc Willotd ||2004 hereafter GRW) in 
which we presented a method of using generalized luminos- 
ity functions (GLFs) to predict the number of radio galax- 
ies and radio-loud quasars as a function of redshift z, radio 
luminosity, narrow-emission-line luminosity and type of uni- 
fied scheme. The GRW analysis introduced two new param- 
eters: a encoding the well known, and strong, correlation 
between radio and narrow-line luminosity; and /3 encoding 
the scatter about this relation. GRW concluded that the 
simplest unified schemes, in which the probability of a given 
object being viewed as a quasar is independent of q, are 
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strongly ruled out by the data available on radio-selected 
samples. However, their analysis was unable to discriminate 
between the two most popular modifications to the simplest 
schemes, namely: (i) models in which the probability of a 
given object being viewed as quas ar is dependent on a and 
0, as in the 'receding torus' scheme iLawrencellQQltlSimpsonl 
Il998r : and (ii) models in which there are two populations 
of radio sources, only one of which necessarily comprises 
of an adm ixture of unobscured and obscured quasar nuclei 
iWillott et a l. 2000') . Both modified schemes were able to ex- 
plain both the systematic increase in the fraction of quasars 
with a and the systematically higher narrow-emission-line 
luminosities for radio quasars over radio galaxies, matched 
in radio luminosity and redshift. 

Differences in luminosities between samples of radio 
quasars and radio galaxies, matched in radio luminosity 
and redshift, have also b een found at submill imetre and in- 
frared wavelengths (e.g. IWillott et al]l2002ah . These wave- 
bands would seem to be ideal for testing unified schemes. 
The putative opaque dust torus is a vital ingredient of all 
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unified schemes as it is needed to hide the active nucleus 
by absorbing the soft X-ray, UV, optical and near-infrared 
radiation. Much of this absorbed energy will be re-radiated 
by dust grains at infrared wavelengths, and for high-redshift 
objects this may be observed in the submillimetre waveband. 
This dust is expected to be optically thin (at least at longer 
wavel engths) and its thermal emission is r adiated isotropi- 
cally (Efstathiou & Rowan- Robi nsonlll994) . 

■Willott et al . (2002a) found that radio- loud quasars 
were a factor of approximately 2 — 5 times brighter than 
radio galaxies from SCUBA data at 850 ^m, where the 
large uncertainty arises because of the difficulties involved 
with correcting the submillimetre flux densities for syn- 
chrotron contamination. However, earlier studies of far- 
infrared differences between quasars and radio galaxies led 
to conflicting results. Broad-lined objects were found to 
be more luminous than narrow-lined objects in the far- 
infrared from fo ur pairs of objects in the redshift rang e 
0.25 < z < 0.57 ivan Bemmel Barthel fc de Graauwll200(Tl . 
whereas iMeisenheimer et al.l ll200ll) found no difference in 
the mid- and far-infrared flux densities of ten pairs of radio 
galaxies and radio quasars with z ^2. Such results have led 
to discussions of the viability of unifled schemes for radio- 
selected quasars and galaxies. 

In this paper, submillimetre data on radio-selected sam- 
ples will be incorporated into the GLF approach of GRW. 
The availability of submillimetre data for complete sam- 
ples of radio sources is much more limited than is the case 
for narrow emission-line data, and non detections are far 
more common. However, the available submillimetre data 
has been collected in a way that yields quantifiable biases. 
We have four main aims. First, to discriminate robustly be- 
tween the competing unified models that were considered in 
GRW. Second, to quantify the correlation between a and 
submillimetre (850 /.im) luminosity Lgso, leading to an es- 
timate of the correlation between (the unobscured) optical 
quasar luminosity Mb (or, equivalently, bolometric quasar 
luminosity Lboi) of radio sources and Lgso. Third, to obtain 
increased predictive power for the GLFs due to the addition 
of submillimetre data. Fourth, to understand the discrepan- 
cies between previous studies of differences between radio 
galaxies and quasars in the far-infrared and submillimetre. 

This paper is organized as follows. In Sec. |21 the data 
to be used are described. In Sec. |3 the correlations be- 
tween submillimetre luminosity and the principal compo- 
nents found in GRW are investigated. In Sec. 21 new GLFs 
based on this correlation are described and the results and 
analysis are presented. We make comparisons between the 
predictions of our GLFs and previous studies in Sec. |5] In 
Sec. |S| we derive a scaling relation between far-infrared lu- 
minosity and the bolometric luminosity of the quasar nu- 
clei of radio sources. We assume throughout that Ho = 
70 kms~^ Mpc"\ r^M = 0.3 and = 0.7 



2 DATA 

A sample of eighteen high-redshift (1.3 < z < 3.5) radio 
quasars with existing 850-/im SCUBA data has been de- 
fined. These are drawn from a suite o f low-frequency-selected 
redsh ift surveys: the 3CRR sample ([La ing. Ril ev fc Longaiil 
Il983l) . the 3CR sample JSpinrad et alJ Il985l) : the 6CE 



sample jRawlines. Bales fc Lacl 1200 ill: parts I and II of 
the 7CRS sample dWillott et alJ l2002bl) : and the TOOT 
(TEXOX-1000) Survev dHifi fc Rawlings! l2003ft . The first 
four of these surveys have complete optical identifications 
and essentialy complete spectroscopic redshifts (or, in the 
case of the 7CRS, reliable photometric redshifts for the few 
objects without secure spectroscopic values). The subsec- 
tions of the TOOT survey used in this work have complete 
spectroscopy but sign ificantly lower redshift completeness 
jHiU fc Rawlingj|2003h : we are confident, however, that few, 
if aiiy, quasars with R 5C 23.5 have been overlooked. 

IWillott ct al. ( 2002a) presented SCUBA data for eleven 
quasars in the redshift range 1.37 < z < 2.0: two 
quasars from the 6CE sample with 2.0 ^ Sisi 3.93 Jy 
(6C 0955-^3844 and 6C 1045+3513); seven quasars from 
the 3CRR sample, within the Right Ascension (RA) range 
04'' < RA < 22''; and two quasars (3C 298 and 3C 280.1) 
from the 3CR catalogue, which were not included in the 
3CRR for reasons of low declination and (slightly) low ra- 
dio flux-density respectively. The sky area for the combined 
3CRR/3CR subsample was estimated by taking of the 
sky area that the 3CR survey would have with the reduced 
RA range, since all of the 3CRR objects are included in the 
3CR survey, but one 3CR object (3C 418) is in the correct 
redshift and RA range but could not be considered here be- 
cau se it has no SCUBA d ata. 

iRawlings et alJ (|20o3) presented submillimetre data for 
eight quasars with 2.5 ^ z 3.5: flve radio quasars se- 
lected from the TOOTOO and TOOT08 sub-regions of the 
TOOT with R s; 23.5 and 0.1 < Sisi < 0.2 Jy, where 
5*151 is the 151-MHz flux density; and three radio quasars 
from parts I and II of the 7CRS sample with 5*151 > 0.5 
Jy. We have halved the effective sky area of the 7CRS sur- 
vey to account for the fact that only two out of four 7CRS 
quasars in the redshift range are used; this assumes that the 
excluded objects, synchrotron-dominated quasars (5C 6.288 
and 5C 6.291; one of which was observed by Rawlings et al.), 
are unbiased with respect to their (isotropic) thermal dust 
luminosities. 

The [OIII] luminosities of the TOOT objects were es- 
timated from the Lyman-a line luminosities (which are 
dominated by emission from t he narrow-line region) from 
iRawlings. Hill fc WiUotd (|2004) and converted t o L[oni] us- 
ing the line ratios presented in lMcCarthvl (|l_993') . [There are 
uncertainties in this procedure because reddening can have 
a large impact on the Lya / [OIII] ratio.] 

For the two 7CRS objects, the [OIII] line luminosity 
was estimated by assuming a rest-frame equivalent width of 
30 A for this line. For the 6 CE and 3CR o bjects, the line 
luminosities are taken from IWiUott et alJ J2002al) . except 
3C437 which had its [OIII] luminosity estimated using the 
[OIII]/[OII] relation from GRW, and 6C 0919-^3806 which 
now c orrectly employs the CII]/[OIII] ratio of McCartl^ 
Jl993h . 

Table shows the basic data for the sample, includ- 
ing the a and (3 values derived following the prescription of 
GRW. Submillimetre luminosity is used as a measure of the 
radiation from dust, which may be heated by star-formation 
and/or the central quasar nucleus. However, since the 6CE 
and 3CR samples include high-flux-density radio sources, 
it is necessary to attempt to correct for contamination of 
the submillimetre radiation by the synchrotron radiation 
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Table 1. Submillimetre observations and basic data for the sample of eighteen radio quasars froml Rawlings et al.l l2004^ and lWillott et all 
{2002^. The values for the [OIII] line strength have been calculated by various methods (see Sec. |2J. Submillimetre detections at the 
^ 2a level are given in bold type. SC is the assumed fractional synchrotron contamination at 850 /^m. D is the projected linear size of 
the radio source and a f denotes that is has been measured between the core and an ext ended component on j ust one side of the core. 
The synchrotron contaminations, expressed as fractions of observed 5850, were taken from lRawlings et alJ i2004h . The parameters a and 
/3 have been derived using the method of GRW. 



responsible for radio emission. The fractional corrections 
for synchrotron contamination, included in Table [J for the 
6CE a nd 3 CR sub-sample s , were taken from .Rawlings et alJ 
J2004fl and IWillott et alJ ^0^^. Throughout this paper, 
submillimetre flux densities are converted to Lgso by, first, 
correcting for synchrotron contamination and then, second, 
by assuming a grey body emission spectrum with a dust 
temperature Td = 40 K and emissivity index /3d ~ 2 
JPriddev fc McMahonlliool . 

Fig. shows the objects in the sample as a function 
of 1/151 and z, in comparison to the full 3CRR/6CE/7CRS 
dataset used in GRW. In this paper we are clearly focusing 
on smaller samples at higher redshifts, but we still span a 
wide range of L151 because of the inclusion of TOOT data. 
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eigenvectors e^/, e^/ and e^i from the PCA analysis, such 

that e^/ = 0.593 logio(Ll5l)norm + 0.571 logio(i[OIII])norm + 

0.569 logiQ(L85o)norm etc. The components are normalized in the 
same manner as GRW and primes are used to distinguish these 
principal components from those found in their previous, two- 
variable, PCA analysis. 



3 THE L850 a CORRELATION 

GRW encoded the I/[oiii] — ^151 correlation for radio sources 
using the parameters a and /3, derived from a principal com- 
ponents analysis. In this section the correlations between 
submillimetre (850-/im) luminosity 1/850, i'i5i, i[oiii] s^nd a 
are investigated. The aim is to find the best way of incorpo- 
rating submillimetre data into the GLFs derived by GRW. 

3.1 New PCA Analysis 

A new principal components analysis (PCA) was carried out 
to determine the nature of the correlation between L151, 
i[oiii] S'lid 1/850- This analysis used the data of Sec.|21 Al- 
though the overall fraction of submillimetre detections is low 
amongst the fainter radio sa mples, such objects h ave been 
detected in a statistical sense iRawlings et al.ll2004^ . There- 
fore, we replaced any synchrotron-corrected values of S85o 
below the la detection limit by this limit for the purposes of 



calculating L850 . The mean values of logj^Q L151 , logj^Q I/[oiii] 
and logj^Q 1/850 are 27.73, 36.65 and 22.68 respectively, with 
associated standard deviations 0.60, 0.55 and 0.38. The re- 
sults of the PCA are given in Table |5| The dominant cor- 
relation (accounting for ~ | of the scatter) is a positive 
correlation between 1/850, i[oiii] a-nd I/151. To test if this 
Icr-upper-limit approach yielded serious problems, we re- 
peated the analysis with ten-times lower values of logj^Q Lsso 
for the non-detections, which yielded only small changes in 
the derived values. Note, however, that this would increase 
the dispersion in logjQ 1/850, so that the interpretation of the 
normalized quantities would be different. 

3.2 Survival analysis 

The new PCA (Sec. 13.111 has yielded the basic result that 
^850 is positively correlated with Z/151 and L[oni]- As the a 
parameter from GRW encodes the L151 and L[oni] correla- 
tion, it is expected that Lgso will be correlated with a. In 
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Figure 1. The 151 MHz luminosity Lisi versus redshift z plane for the 3CRR, 6CE and 7CRS complete samples with the symbols 
identifying radio quasars and radio galaxies from different samples: 3CRR quasars (filled squares); 3CRR radio galaxies (open squares); 
6CE quasars (filled triangles); 6CE radio galaxies (open triangles); 7CRS quasars (filled circles); and 7CRS radio galaxies (open circles). 
Weak quasars (see GRW) are identified by star symbols. Also shown are the five TOOT quasars (inverted triangles), and the two SCR 
quasars (plus signs) described in Sec. |^ Circles surrounding symbols show objects observed by SCUBA, and additional larger circles 
show objects detected by SCUBA at the > 2a level. The dotted line shows the flux-density lower limit for the 3CRR sample, the dashed 
lines show the upper and lower limits for 6CE, the solid line shows the lower limit for 7CRS and the dot-dashed lines show the upper 
and lower limits for the TOOT sample (all assuming a radio spectral index of 0.8). 



this section, a survival analysis is used to quantify the cor- 
relation between log^g Lgso and the logjQ L151 , log^„ Lfomj 
combination through a, now taking into account the 2a up- 
per limits which are important because nearly half of the 
objects are not detected at the 2a level at 850/^m. A cor- 
relation and reg ression analysis was carried out using the 
AsURV software Jlsobe. Feigelson fc Nelsorj[l98^ . 

Table |3 shows the correlation coefficients using three 
methods; Cox's Proportional Hazard model, the generalized 
Kendall's rank correlation test for censored data sets using 
the BHK test, both of which measure the degree of associa- 
tion of two variables, and the Buckley-James method which 



is the usual least squares regression adapted for censored 
data. Fig. |5| shows the correlations of Lsso with L151 and 
L[oiii], and also the correlation of Lgso with a. 



This battery of statistical tests suggets that there are 
real, positive correlations between the variables considered. 
However, the best-fit lines are strongly dependent on how 
the upper limits are treated. In the next section we will ex- 
tend the GLFs from GRW to encode the Lgso — a correlation 
and to treat upper limits in a way more appropriate to the 
problem under consideration. 
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Figure 2. The correlations of logj^Q Lsso with logj^Q L151, logjQ i[oiii] a^nd a for the sample of eighteen quasars defined in Table ITl 
The 2(7 upper limits are shown for the non-detections (symbols as in Fig.0. The best-fitting (dotted) lines are shown from an analysis 
using the Buckley-James method. 
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Table 3. The correlation of logj^Q Lgso with logj^QLisi, logj^g ^[OIII] ^^'^ ^ using Cox's Proportional Haz ard model, the generalized 
Kend all's rank correlation test (BHK test) and Buckley-James method for regression with the Asuav program llsobe. Feigelson &: NelsonI 
[Ss^. is the global x^i and Pc is the probability that a correlation is not present for Cox's Proportional Hazard model, r is Kendall's 
T, Z is the Z-value and Pk is the probability that a correlation is not present for the BHK test. SD is the standard deviation of the 
gradient and SDR is the standard deviation of regression for the Buckley— James method. 



4 GENERALIZED LUMINOSITY FUNCTIONS 
4.1 Definition 

The logical extension to the GLF pq(Li5i, L[oni], 2) = 
pQ{a, 13, z)\J\, where J is the appropriate Jacobian defined 
by GRW, is to incorporate submillimetre data and define 
a GLF cf) which predicts the number of radio quasars per 
unit comoving volume (in units of Mpc'^) per unit base-ten 



logarithm of I/151, L 



[OIII] 



and 1/850 at a given redshift 



The PCA (Sec. 13.111 and survival analyses (Sec. 13.211 have 
demonstrated that the Z/gso-Q correlation can be encoded 
as a linear relation logj^Q Ls5o{q) = Ga + C, where G is the 
gradient and C is the intercept, plus scatter about this re- 
lation. We encode the scatter using a Gaussian distribution 
fiinction about Ls5o{a), with a free parameter aL, so that 



<^(Li5l, L[oiII], 1/850, 



dlogio Lisidlogio L[oni]dlogi(, LgsodV" 
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and dV represents the comoving volume element at an epoch 
corresponding to redshift z. 



Given the large proportion of data with only upper lim- 
its at 850 nm, we need to adapt this GLF so that it takes 
non-detections properly into account. We will assume that 
an observer will measure a flux density Sobs at 850 ^m (after 
correction for synchrotron contamination) which is drawn 
from a Gaussian distribution about the actual fiux density 
5*850 with the standard deviation given by the observational 
error cTobs- Thus, an observed value of flux density Sobs could 
be due to a range of real (synchrotron-corrected) fiux den- 
sity S850, each of which corresponds to a lumin osity Laso 
(assuming /3d and Td from IPriddev fc McMahonll2001^ ■ In 
this way, we define a new GLF $ in terms of Sobs and ob- 
servational error CTobs, related to the GLF (j) by 



'l>(I/151,I/[OIII] z) (3) 

= J </'(-f'151,-f'[OIII],-^'850,2)<7(Sobs,S850,CTobs)iJ''idS850, 

where 



5(Sobs, S850, CTobs) ~ Prob(Sobs|S850, CTobs) 



(4) 



1 / [Sobs — Sg5o] 

exp ' 



CTobs VStT 



2<^o\s 



and |J'| is the new Jacobian needed to convert from S85o to 
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4.2 Normalization corrections 

One of the known small problems with the GLFs defined by 
GRW is that they predict too few objects at low redshift, 
and correspondingly too many objects at high redshifts in 
comparison to the 3CRR, 6CE and 7CRS data. This was a 
consequence of trying to fit the 6C and 7C source counts (in 
addition to the 3CRR, 6CE and 7CRS redshift surveys) and 
may reflect imperfections in the evolutionary model used 
for the high-a population (GRW). Since the redshift range 
of our sample is 1.37 < z < 3.5, and we have only a small 
sample of objects, we expect that such imperfections may be 
important. Table|3]shows how we have attempted to correct 
for this effect by scaling the sky areas of each subsample by 
the ratio of the total number of objects in each subsample 
to the nearest whole number of objects (quasars and radio 
galaxies) predicted by the unified schemes for each of the 
models (2RS, 2SS and 2RR; see Table 4 of GRW). Note 
that we do not yet know the number of radio galaxies in the 
TOOT surveys in this redshift range, so we have estimated 
the correction factor by using the average quasar fraction, 
from the three models, to predict the number of TOOT radio 
galaxies from the number of TOOT quasars. 

4.3 Results 

The GLF of Eqn. Q was combined with each of the three 
models (2RS, 2SS and 2RR) defined in Table 4 of GRW 
using the best-fitting unified scheme parameters found in 
their Table 5. A maximum-likelihood analysis was performed 
in order to find best- fitting values of ctl, G and C. The 
likelihood function rj was defined as 

18 

V = —2 ln[<I>(Li5i . , -L[oiii]; , 5'obsi , (Tobsi , Zj)] 

1=0 

+ 2 JJJJ[<j){Li5i,L[oiii],Ls50,z)Q{Li5i,z) 
dV 

X dzd(logi(, Li5i) d(logi(, I/[oiii]) d(logi(, L85o)],(5) 

where is the sky area available for each sample of radio 
sources, modified by the factors derived from Table |1] 

The best- fitting values of the three free parameters a^, 
G and C were found to be almost identical for each of the 
three models and are presented in Table |S| The steeper gra- 
dient and different normalization of this line compared to 
that found by the survival analysis can be explained by our 
new method of treating upper limits which means that non 
detections can be drawn from a wide range of underlying 
submillimetre luminosities. 

A probability distribution was found for the gradient G 
of the Lgso versus a correlation by marginalising the like- 
lihood (the function t; in Eqn.|SJ over C and ctl, and nor- 
malizing. Fig. 13 shows that the probability that there is a 
positive correlation {S > 0) between Lsso and a is ~ 99 
per cent. Note that this is a different conclusion than one 
would draw from the 'formal' error from Table |3 because 
marginalization allows correctly for the complicated shape 
of the probability distributions. Note, however, that there is 
still quite a large error on the gradient G of the correlation. 

The relative probability of the models with a receding 



0.12 r 



0.10 - 
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Figure 3. Marginalized probability distributions for the gradi- 
ent G of the logj^Q Z/850 — correlation derived from the 2RS 
GLF (similar probability distributions can be inferred using mod- 
els 2RR and 2SS). Different datasets have been used: the radio 
quasars from TablelTl fdot-dashed line); dataset A (dotted line); B 
(solid line); and C (dashed line). Datasets A, B and C are defined 
in Sec. O . 

torus in the high-Q population (models 2RS and 2RR) with 
respect to the GLF with a non-luminosity-dependent unified 
scheme in both populations (model 2SS) is ~ 100. The main 
reason for this is that the adoption of a receding torus means 
that the value of the transition angle Otrans increases with 
Q, and as all the radio quasars in the sample have high val- 
ues of a (Fig. |5J , Eqn. 10 of GRW implies a corresponding 
boost in the values of pQ [note that the free parameter fix- 
ing the behaviour of Otrans has not been allowed to vary in 
our new likelihood analysis]. However, in complete samples, 
these increased likelihoods of seeing radio quasars goes hand- 
in-hand with decreased likelihoods of seeing radio galaxies 
(see Fig. 8 of GRW), so we defer interpretation of any dif- 
ference in likelihoods until, in Sec. 15. II we consider complete 
samples including both radio quasars and radio galaxies. 



5 COMPARISON WITH PREVIOUS STUDIES 

5.1 Comparison with the Willott et al. sample 

IWillott et alJ l|2002^ found that radio quasars were a factor 
of ~ 5 times brighter than radio galaxies in the submillime- 
tre, based on SCUBA observations of eleven 3CRR, 3CR 
and 6CE quasars (see the lowest rows of Table0 and twelve 
3CRR and 6CE radio galaxies (Table ij. This factor could 
be reduced somewhat by adopting larger estimates of syn- 
chrotron contamination, but not below a factor ~2, and so 
not to the point where radio quasars and radio galaxies have 
comparable dust luminosities. We next make a comparison 
of the ability of the three GLFs to reproduce the distribu- 
tion of submillimetre fiuxes of these data. A one-dimensional 
optimization algo rithm (the Golden Section Search method; 
iPress et al.lll992^ was used to find the best-fitting value of 
itl , using the data (on both radio qu asars and radio galaxies) 
presented bv lWillott et alJ (l2002af) . As these radio quasars 
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RG 
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RG 
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RG 
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RG 
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RG 


Data 


9 


7 


2 


5 


2 





5 


7 


18 


7 


2RS 


12.38 


8.23 


7.13 


5.18 


2.64 


1.92 


7.88 


21.24 


32.05 


38.04 


2SS 


8.24 


12.23 


7.71 


5.39 


2.58 


1.87 


11.37 


17.75 


31.06 


39.27 


2RR 


12.65 


7.97 


7.41 


5.37 


2.58 


1.88 


8.26 


20.86 


33.12 


37.19 


2RS 


9.43 


6.27 


4.53 


2.31 


0.78 


1.28 


4.35 


11.72 


19.07 


21.99 


2SS 


6.28 


9.32 


4.72 


2.26 


0.81 


1.25 


6.27 


9.79 


18.37 


22.85 


2RR 


9.64 


6.07 


4.70 


2.26 


0.81 


1.25 


4.56 


11.51 


19.70 


21.49 



Table 4. Table of the number of radio quasars (Q) and radio galaxies (RG) from each survey in the sample (upper panel) and the 
predicted numbers of quasars and radio galaxies in each subsample before and after (middle and lower panels respectively) a factor was 
applied to the sky areas to correct for the imperfections in the GLFs at high redshifts (Sec. see GRW for a full description of the 

three models (2RS, 2SS, 2RR) which represent the three different unified schemes. Note that the number of radio galaxies in the relevant 
redshift range is not yet established for the TOOT survey. 



Model 




G 


C 




dct(VV»7min) 


-P2RS 


2RS 


U.dDO_o.ioi 


6.036 ±5.180 


22.136 ± 0.493 


651.79 


1534.28 


1.000 


2SS 


U.ciDO_(, -^p^ 


6.036 ±5.180 


22.136 ± 0.493 


661.37 


1534.27 


0.009 


2RR 


'^■■^"0-0. 101 


6.036 ±5.180 


22.136 ± 0.493 


652.04 


1534.28 


0.886 



Table 5. Results of the likelihood analysis of Sec. 14.31 which is based on the sample of radio quasars in Table The best-fitting 
values of the parameters , G and C, the minimum value of the likelihood function ?7min corresponding to these best-fitting values and 
det(VV»;niin) for the three unified schemes. 'Formal' error bars on the parameters are obtained from the VVj7 matrix in the normal 
manner (see GRW), and the probability P2RS of each of the models with respect to model 2RS is calculated using Equation 12 of GRW. 



Model g-L ^min det(VVr?min) P2RS 

2RS QASl+'^^Xlo 876.88 isl 1.000 

2SS 0.48lto.}4o 884.57 45.2 0.028 

2RR 0.481+n?^I 877.45 45.2 0.752 



Table 7. Results of the likelihood analysis of Sec. I5.ll which is 
based on t he radio quasars and radio galaxies in the sample of 
IWillott et a l. 1 2002a). The best-fitting parameter (T^, the mini- 
mum value of the likelihood function )7niin corresponding to this 
value and det( VVr/jjiin) for the three unified schemes. 'Formal' 
error bars on the parameters are obtained from the VVr; matrix 
in the normal manner (see GRW), and the probability P2RS of 
each of the models with respect to model 2RS is calculated using 
Equation 12 of GRW. 



a significant diflerence in the Lsso values of radio quasars 
and radio galaxies. For this sample, we find that models 
with a receding torus in the high-a population (2RS and 
2RR) are strongly favoured (by a factor of ~ 40) over a 
model with a torus which has a constant opening angle in the 
high-a population. This constrasts with the very marginal 
(factor ~ 2) difference in likelihoods between such models 
in GRW. Since, essentially, the L151 and L[oni] data con- 
sidered are a subset of the data considered by GRW, and 
(in contrast to the likelihood analysis of Sec. 14.31 both ra- 
dio quasars and radio galaxies are included, the reason that 
we are reaching such a different conclusion is that this cur- 
rent study focuses on luminosity differences between radio 
quasars and radio galaxies in just the high-a population. 
The use of a relatively narrow range of redshifts, so that 
the confusions associated with mixing low- and high-a pop- 
ulations and with wide ranges of redshift is an additional 
factor, as is the additional cons traint of the submillimetre 
properties. IWillott et al.l ll2002al) attributed the difference in 
submillimetre flux-densities between radio quasars and radio 
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Figure 4. The 850/xm luminosity Lgso against a for the sample 
of eleven quasars (last 11 r ows of Tablepi and tw elve radio galax- 
ies (Table|HJ considered bv lWillott et alj ^2002a^ with symbols as 
in Fig. 2cr upper limits are shown for the non-detections. The 
best-fitting line from the optimization procedure is shown (solid 
line; see Sec l5.1l . 



and radio galaxies now span a much smaller range in a, 
compared to the full sample used to constrain the gradient 
G and offset C (Sec. 14.31 . we initially kept the best-fitting 
values of G and C constant during this optimization. 

We present the results in Table |7] Fig. 01 shows this 
sample with the best-fitting gradient determined by the op- 
timization procedure in Sec. 14.31 Note that the addition of 
the radio galaxies leads to a small increase in the derived 
scatter aL in the relationship between Lsso and a, but that 
the scatter is still insufficient to hide what appears to be 



8 J. A. Grimes, S. Rawlings & C.J.Willott 



Name 


Class 


Redshift 


D/ 


logic) [-f-isi/ 


logio[-^[Oni]/ 


Ss50/ 


5450/ 


SO 


a 


P 








z 


kpc 


WHz-isr-i] 


W ] 


mJy 


mJy 








4C 13.66 


G 


1, 


.450 


50.7 


28.09 


36.63 


3.53 ± 0.96 


-16.20 ± 18.20 


0.00 


0.095 


0.010 


3C 329 


G 


1, 


.781 


94.6 


28.42 


36.38 


0.83 ± 1.00 


-2.10 ± 18.30 


0.00 


0.099 


0.031 


3C 241 


G 


1, 


.617 


7.7 


28.11 


36.75 


1.81 ± 0.94 


14.80 ± 12.60 


0.00 


0.101 


0.006 


3C 294 


G 


1, 


.786 


126.7 


28.39 


36.97 


0.19 ± 0.78 


5.40 ± 13.40 


0.00 


0.120 


0.008 


3C 322 


G 


1, 


.681 


279.4 


28.21 


35.97 


-0.05 ± 1.06 


-37.50 ± 16.00 


0.00 


0.076 


0.038 


3C 437 


G 


1, 


.480 


291.0 


28.11 


36.37 


-1.18 ± 0.98 


2.90 ± 17.30 


0.00 


0.086 


0.020 


3C 470 


G 


1, 


.653 


203.2 


28.15 


35.93 


5.64 ± 1.08 


57.80 ± 32.90 


0.00 


0.072 


0.037 


6C 0820+3642 


G 


1, 


.860 


193.8 


27.55 


36.01 


2.07 ± 0.96 


13.60 ± 18.00 


0.00 


0.052 


0.011 


6C 0901+3551 


G 


1, 


.904 


21.9 


27.46 


36.13 


-1.83 ± 1.15 


-19.30 ± 8.20 


0.00 


0.053 


0.003 


6C 0905+3955 


G 


1, 


.882 


934.5 


27.71 


35.77 


3.62 ± 0.89 


31.20 ± 16.20 


0.00 


0.049 


0.026 


6C 0919+3806 


G 


1, 


.650 


88.1 


27.53 


35.93 


-0.88 ± 1.05 


10.50 ± 10.10 


0.00 


0.048 


0.014 


6C 1204+3708 


G 


1, 


.779 


435.2 


27.61 


36.45 


0.16 ± 1.25 


45.10 ± 26.60 


0.00 


0.070 


-0.003 



Table 6. Submillimetre observations and basic data for the sample of twelve radio galaxies from lWillott et alj (2002^. The values for 
the [OIII] line strength have been calculated by various methods (see Sec.|5J. Submillimetre detections at the ^ 2cr level are given in 
bold type. SC is the assumed fractional synchrotron contamination at 850 /xm. D is the projected linear size and a f denotes that it is 
measured between the core and an extended component on just one side of the core, contaminations, expressed as fractions of observed 
S'gso, were taken from iRawlings et al.l l2004h : o and /3 are the parameters defined by GRW. 



galaxies to the dust being heated by quasars (rather than by 
starbursts) which would cause emission-line luminosities to 
be correlated with submillimetre luminosities if both reflect 
the luminosity of the central source. This result is consistent 
with a receding torus because objects with more luminous 
emission lines are more likely to be viewed as quasars, be- 
cause of the larger torus opening angles, and therefore to 
have higher I/[oiii], a and I/gso- 

The weighted mean of the synchrotron- corrected sub- 
millime tre fluxes for the radio quasars in the lWillott et afl 
sample is 4.72 + 0.36 mJy, and for the radio galaxies 
it is 1.25±0.29 mJy; the median is 4.47 mJy for quasars and 
0.51 mJy for radio galaxies. We have simulated distributions 
of Sgso from the average of 30 simulations using each of the 
three GLFs. For model 2SS, the shift is ^ mJy with the 
median Ssso being 3 mJy for both radio quasars and radio 
galaxies. This is as expected because, without a receding 
torus, experiments that focus only on the high-a population 
would find no difference in the a and hence Sgso values of ra- 
dio quasars and radio galaxies. We find that the shift in 5*850 
between the median quasar and the median radio galaxy is 
1.0 mJy for models 2RR and 2RS (the median quasar has a 
fiux density of 3.0 mJy and the median radio galaxy has a 
fiux density of 2.0 mJy), so that quasars are predicted to be 
1.5-times brighter in the submillimetre than radio galaxies. 
Encouragingly, then, shifts in the direction seen in the data 
are predicted, but not by a large enough factor. It is easy 
to see how this factor has arisen. The median a is 0.113 
for quasars and 0.074 for radio galaxies due a larger me- 
dian emission line luminosity of logj^Q(L[oiii]/W) = 37.00 
for quasars than the corresponding value of 36.25 for ra- 
dio galaxies, with relatively well-matched radio luminosi- 
ties [median logio(Li5i/(WHz~^sr~'^)) for quasars of 28.19 
compared to 28.10 for radio galaxies]. Using the best-fitting 
^850 — Of relation, this difference in a leads to median 
logio(L85o/WHz~^sr"^) of 22.82 for quasars and 22.58 for 
radio galaxies - a difference of 0.24 or a factor of « 1.7. 

It follows that the optimized gradient for the strength of 
the L850— Q correlation may not be hi gh enough to reproduc e 
the submillimetre difference seen bv lWillott et alJ (|20023). 
However, this gradient was derived based on quasar data 



alone, and we have seen from the error bars given in Table 
|5] and the width of the marginalized probability distribution 
in Fig. |3] that there is a large uncertainty in these values. 
We therefore repeated the optimization procedure with three 
free parameters [cjl , G and C) and for three datasets defined 
as follows: (A) the eighteen quasars from Table Q and the 
twelve radio galaxies from Table |S| (B) the eleven quasar s 
and the twelve radio galaxies from IWillott et al.l ll2002ar) ; 
and (C) the same dataset as (A), but excluding one radio 
quasar, 6C 1045+3513, the most prominent outlier in Fig|^ 

In Table |S] we present the results of the optimization 
procedure, and the lines using these best-fitting gradients 
and offsets are plotted in Fig. |S] We find that models 2RS, 
2RR and 2SS all give very similar optimized parameters for 
a given dataset (A, B or C), and that for all datasets mod- 
els with a receding torus (either 2RR or 2RS) in the high-a 
population are favoured strongly over the 2SS model. Formal 
errors could not be found for the gradient G and the offset 
C for dataset B because the range in a for this dataset is so 
small that the likelihood surfaces are not well described by a 
multi-dimensional-Gaussian approximation. To find an error 
on G (given in the brackets in Table|5J, we fitted a Gaussian 
to the marginalized probability distribution of the gradient 
(Fig. 1^ ; the probability distributions for G given the 2RR 
model and the various datasets are very similar. The error 
on C, given in the brackets in Table |H1 is a rough estimate 
based on the expectation that it is greater than that for 
dataset A, because of the small range in a, but not as great 
as that for dataset C, which arises because of the double- 
peaked shape of the probability distribution which is more 
pronounced than that for dataset A (Fig. These pecu- 
liar double-peaked shapes arise because of the large fraction 
of upper limits in the datasets. Note, however, that for all 
datasets, a gradient of zero is strongly disfavoured because 
the upper limits are concentrated at low values of a. Note 
also, that the secondary peaks at 5* ~ 20 are strongly dis- 
favoured by the statisti cal detection of the TOOT quasars 
in the submillimetre bv lRawlings et al] j|2004h . so ignoring 
these, the peaks in the marginalized probability distribu- 
tions all cluster in the range 6^; 10. 

We see in Fig. |3 that the large uncertainty in the gra- 
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Dataset 


0"L 


G 


C 


2SS: P2RS 


2RR: Paas 


A 
B 
C 


475+"-^'*« 

'^•^''^-0.114 

r. 000 + 0.169 

u.oao_Q j^j^g 

U.4dS_(, 


6.64 ± 7.90 
9.91 ± (7.0) 
9.41 ± 14.66 


21.89 ± 0.72 
21.56 ± (1.0) 
21.61 ± 1.39 


0.018 
0.034 
0.015 


0.989 
0.898 
1.022 



Table 8. The best-fitting parameters tr^, G and C for datasets A, B and C (see Sec. 15. ll for a definition of these datasets). 'Formal' 
error bars on the parameters are obtained from the VVr; matrix in the normal manner (see GRW), and the probability P2RS of each of 
the models with respect to model 2RS is calculated using Equation 12 of GRW. Brackets indicate that the error is an estimate, because 
formal errors could not be found for these parameters (see Sec. l5.ll for a discussion of how these were estimated). 
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Figure 5. Graphical representation of the results of the optimization procedure described in Sec. lOJand Table ISl The lines showing 
the best-fitting values for the Lgso ~ correlation for datasets A, B and C, are shown alongside the relevant data (symbols as in Fig.lTl. 



dient G is also compounded by the dramatic influence of 
just one outlying point (6C 1045-1-3513). It is interesting 
that the peak of the probability distribution for G shifts to 
significantly higher values for datasets B and C, in which 
this prominent outlier is excluded. It is probable that the 
submillimetre emiss ion for this object i s not dominated by 
quasar-heated dust JWillott et al .l2002al) . but dominated in- 
stead by starburst-heated dust. The small size of its radio 
jets indicate that it is a young radio source in which in- 
tense star-formation activity may be synchronized with the 
jet-triggering event. This hints that the correlation may be 
strongest for more typical radio sources in which quasar- 
heated dust may dominate even at the longer wavelengths 
probed in the submillimetre. 

Simulations of the Sgso — iisi plane using the opti- 
mized parameters from dataset C were performed. Once 
more, as expected, the 2SS scheme gives identical values 



of Sgso (~ 3.40 Jy) for radio quasars and radio galaxies. For 
the other schemes, larger shifts are confirmed: an average ob- 
served fiux Ss50 of 4.65 Jy for radio quasars and only 2.09 Jy 
for radio galaxies for 2RS; and 4.39 Jy for radio quasars and 
2.20 Jy for radio galaxies for 2RR. Moreover, Fig. |3 shows 
that it is perfectly plausible that the true gradient is large 
enough that the difference in submillimetre luminosities of 
radio quasars and radio galaxies could easily be by a fac- 
tor ^ 2. Accounting for residual unce rtainties in the correc - 
tions for synchrotron contamination JWillott et al.ll2002al) . 
the fac tor ■-^2 — 5 difference in submillimetre luminosities 
seen bv lWiUott eTa l. (2002a) are now in quantitative accord 
with both the 2RS and 2RR models. 
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5.2 Comparison with 'matched-pair' experiments 

We consider in this section whether the GLFs described 
in Sec . |1] c an reproduce the results of [ M^ i^^nhei m^r et all 
i200 l') and 'van Bemmel. Barth el fc de &aauwr(l2000r !^ 

Mciscnhcirncr et al. (2001) observed ten pairs of SCR 
radio galaxies and radio quasars with ISOPHOT on ISO 
at infrared wavelengths between 5 and 180/im. The sources 
were matched in redshift and 178-MHz luminosity Lirg, but 
spanned a large redshift range of 0.05 < z < 2.02. The av- 
erage value of 100-micron flux density 5ioo was 293 mjy 
for quasars and 300 mJy for radio galaxies, leading to their 
conclusion that quasars and radio galaxies cannot be distin- 
guished by their observed mid- and far-infrared properties. 
However, if we use /3d ~ 2 and the val ues of dust temper- 
atures found by fitting to the data in iMeisenheimer et alJ 
i200lfl [if no values were given, a value of 40 K was as- 
sumed], we find that the average (rest-frame) 100- /im lumi- 
nosity I/ioo was 38.69 for quasars and 38.36 for radio galax- 
ies, so that their radio quasars are an average of ~ 2-times 
more luminous than radio galaxies at 100 /im. 

We simulated histo grams of logj^g -^'loo for the 
IMeisenheimer et al.l i200 j) dataset from the average of 30 
simulations using the optimized parameters for dataset C 
fSec. I5.lt . For 2SS, it is 38.57 for radio quasars and 38.39 
for radio galaxies, with a shift of 0.18 because the range of 
2 is sufficiently large that the low- and high-Q populations 
become mixed, as was the case for the experiments of GRW. 
For 2RS, the median Lioo is 38.75 for quasars and 38.27 for 
radio galaxies giving a shift of 0.48. For 2RR, it is 38.75 for 
quasars and 38.21 for radio galaxies, with a shift of 0.54, so 
that for all models, radio quasars are expected to be brighter 
than radio galaxies, in agreement with the data. To try and 
account for the effects of the selection of matched pairs, we 
used the GLFs to pick a value of Lioo (by extrapolating 
from the value of Lgso using our canonical dust spectrum) 
from the actual red shift and Lirg values of eac h quasar and 
radio galaxy in the IMeisenheimer et alJ i200l[) sample. We 
find that in only ~ 60 per cent of the pairs, is the quasar 
brighter than the radio galaxy in Lioo, based on one thou- 
sand simulations of picking the 10 pairs. This shows that 
the matched pairs method is not very effective at discover- 
ing quite large differences in the properties of radio quasars 
and radio galaxies. ^ 

A similar 'matched pair' experiment was per- 
formed at l ower re dshift (0.25 < z < 0.55) by 
Ivan Bemmel. Barthel fc de Graauwl i2000ll . They obtained 
ISOPHOT data for four pairs of 3CR radio galaxies and 
radio quasars, matched in 178 MHz luminosity and red- 
shift. They concluded that broad-lined objects were brighter 
in the far-infrared than narrow-lined objects, since the av- 
erage 90-/im luminosity Lgo was 37.90 for radio quasars 
and ^ 37.76 for radio galaxies. Simulating their experi- 
ment using our GLFs we could reproduce this measured 
difference between radio quasars and radio galaxies for 
any of the u nified sc hemes. We note that the sample of 
Ivan Bemmel. Barthel fc de Graauw. (■200Qi'l contains such a 

Note also that in one of the lMeisenheimer et al] <200ll) matched 
pairs, 3C 325 (the brighter far-IR source of its pair) is taken to 
be a radio galaxy, whereas our new optical spectroscopic data 
(Appendix A) suggests it is a reddened quasar. 



small number of objects that it would be very difficult to 
reach significant conclusions from such an experiment. 

In conclusion, existing matched-pair experiments typi- 
cally mix low- and high-a populations so all the GLFs pre- 
dict a luminosity difference between radio quasars and radio 
galaxies at 100 /im. We have seen that using small numbers 
of matched pairs is not a particularly good way of compar- 
ing unified schemes, but that the results from such methods 
are consistent with the other results described in this paper. 



6 SCALING RELATION BETWEEN Lgso AND 
Lbol 

Here, we derive a rough scaling relation between Lgso 
and bolometric quasar luminosity LboI for the high-a ra- 
dio source population, argue that it may be applicable to 
quasars in general, and suggest how it might physically arise. 

The form of the underlying Lgso-I/Boi relation can be 
inferred from the gradient G (Table |H] Fig. Using the 
value of G derived from dataset C (see Sec. 15.11 , we find 
that 

a oc ^-^logjoLgso- (6) 

From Equations 1 and 2 of GRW, we have another ex- 
pression for a 

aoc i— logio Li5i + i— logioL[oiii], (7) 

siV2y'n S2\/2y/n 

in which si = 1.05, S2 = 1.10, and n = 302. 
From Equation 4 of GRW, we find that 

logio ^[oiii] oc 1.045 logio Lisi , (8) 

so that if we assume that LboI oc -L[oiii], then 

T r0.7±0.2 /n\ 

Lgso OC Lg^i , (9) 

where the Icr-error bar has been estimated from the proba- 
bility distribution of the gradient (Fig. |3] ignoring the sec- 
ondary peak) rather than from the value from the covariance 
matrix. Note that adoption of dataset A would be consis- 
tent with Eqn. |51 but would favour gradients towards the 
lower end of the quoted error bar. Note also that values de- 
rived from secondary peaks in the marginalized probability 
distributions of G are inconsistent with the constraints set 
from considering the sample consisting of only radio quasars 
[sec Fig.Oand the statistical detection of TOOT quasars by 
Rawlings ct al. (2004)], and seem likely to be artefacts of the 
small ranges in a probed by submillimetre-detected objects 
in the relevant datasets. The scaling relation of Equation 9 
is in good agreement with the results of Sec. 13.11 in which 
it was found that the first principal component represented 
a linear (roughly equal weight) combination of normalized 
versions of log^gLisi, log^g L[o///] and log^Lgso, meaning 
that the three variables scale-up with each other according 
to power laws with indices given roughly by the ratios of the 
associated data standard deviations, namely ~ 1 : 1 : 0.6 
(assuming the real values of logj^Q Lgso lie close to their l-cr 
upper limits). 

The derived Lgso-Leoi scaling relation is simi- 
lar to that inferred from the correlation between 
Lgso and the optical luminosity of radio-quiet quasars: 
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IWillott. Rawlings fc Grimes! J2003l) found a factor ~ 1.4 re- 
duction in Z/850 per unit drop in absolute blue magnitude 
Mb, which corresponds to Lgso oc I/boi if we assume that 
LboI scales linearly with the observed optical luminosity. 
This suggests that the scaling relation derived here is not 
peculiar to radio-selected objects, but is a generic feature of 
AGN. 

There is one obvious way in which this scaling might 
arise. Consider models in which the observed submillime- 
tre luminosity of powerful AGN comes from the dust in 
the obscuring torus as a r esult of d irect quasar heating. 
As suggested previously bv lSimpson fc RawliiiES (2000), a 
receding-torus-like scheme means that the solid angle sub- 
tended by the torus scale with radius r and hence LboI as 
h/r oc -f/goi"' for constant torus half-height h. If the fraction 
of the radiation re-processed by the dust torus is simply 
proportional to this solid angle, then the implied scaling re- 
lation is 1/850 oc I/Boi- This is within the errors of the value 
observed (Eqn.l^J. Of course, more realistic receding torus 
models might feature dis tributions in h and/or systematic 
variations of h with LboI (|simESo3|2003), and hence predict 
slightly different scaling relations 



7 CONCLUSIONS 

In this paper, the generalized-luminosity-function (GLF) ap- 
proach of GRW has been extended to predict the number 
density of low-radio-frequency-selected quasars and galax- 
ies as a function of redshift, radio luminosity, emission- 
line luminosity and 850-/im luminosity, or observed 850-/im 
flux density and sensitivity, using a new method which en- 
ables us to use data with upper flux-density limits. The 
new free parameters in these extended GLFs were fixed by 
comparison with SCUBA (850- /im) data on high-luminosity, 
high-redshift radio sources. Results from our GLF approach 
were found to reproduce previous comparitive studies of the 
submillimetre and far- infrared luminosities of radio quasars 
and radio galaxies IWillott ct al. 2002a; Meiscrihcimer et alJ 
l200ll:lvan Bemmel. Barthel fc de Graauwll20o'^ 

We find robust evidence that the submillimetre luminos- 
ity of low-frequency-selected radio sources is correlated with 
the GRW a parameter, and hence the bolometric quasar lu- 
minosity I/Boi via a scaling relation Lgso oc L'^l^'^'^- This 
corresponds to a factor of ~ 1.9 ± 0.3 reduction in Lgso per 
unit drop in absolute blue magnitude Mb which is similar 
to the correlation tentatively identified for the radio -quiet 
quasar population dWillott. Rawlings fc Grime jl20o3) . The 
fact that the gradient of the scaling relation is less steep than 
a proportionality may give clues as to how Lsso arises, and is 
consistent with models in which quasar heating dominates, 
and in which dust in the receding torus intercepts and re- 
processes a lower fraction of the bolometric luminosity into 
the mid- and far-infrared bands. 

We also find robust evidence that a receding-torus-like 
scheme is strongly favoured (by a factor ~ 40 in probability 
over models without a receding torus) in the high-a popu- 
lation of a two-population scheme for radio sources. In such 
schemes, radio quasars are brighter in narrow emission lines 
and submillimetre luminosity (than radio galaxies matched 
in radio luminosity and redshift) because they are biased to- 



wards the objects with the highest values within the spread 
of bolometric luminosity. 
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APPENDIX A: CORRECTED REDSHIFT AND 
CLASSIFICATION FOR 3C 325 

In this Appendix we present a new optical spectrum of 
3C 325 which has allowed us to correct both the red- 
shift and classification of this object. We targeted this 
object because it is highly unusual amongst high-redshift 
{z > 0.5) radio galaxies in the SCR catalogue at it has 
been detected (using ISO) as a luminous far-infrared source 
iMeisenheimer et al.ll20o'lh . High-z radio sources detected 
in the far-infrared by IRAS and ISO are almost excusively 
quasars, sometimes because the far- infrared emission is from 
a bea med non-thermal component iHes. Barthel_fc_jbsburY 
[^9^V but sometimes because, as detailed "in TwiUot^^^ 
(|2002a|) and this paper, radio quasars seem to be genuinely 
brighter (rest-frame) far-infrared sources than radio galax- 
ies when matched in radio luminosity and redshift. Our 
previous optical spectroscopy of another anomalous far- 
infrared-loud 3CR 'galaxy', 3C 318, found that this object 
had an incorrect redshift and classification in the litera- 
ture, and that the object was actually a z = 1.574 red- 
dened quasar JWillott. Rawlines fc JarvisI |2000D . Prior to 
this work, 3C 325 was classified as an 'N' galaxy with a 
redshift z = 0.86 based on the identification of spectral fea- 
tures with moderat e, possibly broad, M g 2799 and weak 
[Oil] 3727 emission dSoinrad et al.lll985ll . 

3C 325 was observed with the dual-beam ISIS spec- 
trograph on the William Herschel Telescope on 22 April 
2003; the R300B and R316R gratings were used. Condi- 
tions were photometric with good (~ 0.8 arcsec) seeing and 
a 1.5-arcsec slit was used. The telescope was pointed at 
the known optical identification of 3C 325 (RA 15:49:58.22, 
DEC -1-62:41:22.4, J2000.0). The total integration time of 
the observation was 2700s, split into three separate expo- 
sures on each arm. In Figure lAll we show the extracted 
optical spectrum which was spatially unresolved in the two- 
dimensional spectrum. The object has many strong emission 
lines (see Table Al) and therefore an unequivocal redshift 
of z = 1.135; the previously reported redshift is incorrect, 
presumably because the [NelV] line was confused with Mg 
II. 

Although the cores of most of the emission lines are 
narrow there is a definite broad (~ 3000 km s~^) component 
underlying the C III] line, and probably a similar component 
associated with C IV. The Mg II line maybe dominated by 
a broad feature but it is hard to be certain because there 
seems to be a lot of absorption to the blue side of the line. 
The red continuum has a spectral index of « 2.5, so this 
object seems to have all the characteristics of a reddened 
quasar. 

3C 325 is a reddened quasar at z = 1.135, and no t an 
'N' g alaxy at z = 0.86 as previously reported ||Spinrad e t alJ 
[l985|) . It is interesting that so many of the far-infrared- loud 
radio sources at high redshift appear to be in the reddene d 
quasar class, as we have discussed in lWillott et al ]l|20023). 
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Figure Al. Optical spectrum of the reddened quasar 3C 325, with the emission Unes labelled. 
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Table Al. Emission line data for 3C 325 from our optical spectrum. Emission line widths have been deconvolved assuming that the 
source of optical emission was sijatially umcsolved. f denotes that this refers to the core of a line which appears to have broad wings. | 
denotes that values are very uncertain because of absorption. The line fluxes have been measured by integrating along the slit direction 
and no corrections for slit losses have been attempted; estimated errors axe expressed, in brackets, as a percentage of the measured flux. 



